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We consider the possibility of measuring non-equilibrium properties of the current correlation 
functions at high temperatures (and small bias). Through the example of the third cumulant of the 
current {S3) we demonstrate that odd order correlation functions represent non-equilibrium physics 
even at small external bias and high temperatures. We calculate S3 = yieV /T)e^ I for a quasi-one- 
dimensional diffusive constriction. We calculate the scaling function y in two regimes: when the 
scattering processes are purely elastic and when the inelastic electron-electron scattering is strong. 
In both cases we find that y interpolates between two constants. In the low (high) temperature limit 
y is strongly (weakly) enhanced (suppressed) by the electron-electron scattering. 



I. INTRODUCTION 

Since its discovery by Schottky in 1918, noise in electri- 
cal circuits has been thoroughly investigated. Numerous 
studies over the past decade, both experimental^ and the- 
oretical (see Refs.^''^ for reviews), emphasized the quan- 
tum and the mesoscopic aspects of noise, addressing quite 



extensively the issue of non-equilibrium "shot" noise. 
Shot noise has been employed for probing fundamental 
physics of interacting electrons in the FQHE regime^. 
Most of previous work focused on current-current correla- 
tions, i.e., fluctuations derived from second order current 
cumulants^ 



S2(a;,x',w) = i 



d{t ~ t')e"^^*-*") {i{t, x)i{t', x') + i{t', x')i{t, x)), 



(1) 



X 



where /(t, x) denotes the instantaneous current oper- 
ator integrated trough the entire cross-section (at point 
a;); x is the direction of the flow of the d.c. current. Equi- 
librium noise is governed by the Callen-Welton relation, 
also known as Fluctuation Dissipation Theorem (FDT). 
It connects the Ohmic part of the conductance G{uj) with 
the equilibrium current noise (^2 — S2 — ,/=(/)) 



32"^ {uj) — hujG{uj) coth 



huj 
2T 



(2) 



Let us consider a two-terminal constriction subject to 
an applied bias, eV, which is equal to the difference be- 
tween the chemical potentials of the left and the right 
reservoirs eV = — jj,^. The constriction length L 
is much larger than an elastic mean free path I = VfT 
{L ^ I). The value of the d.c. electron current through 
the system is given by Ohm's law, / = GV. In the ab- 
sence of electron-electron collisions the second moment 
of the current fluctuations can be expressed in terms of 
the equilibrium correlation function, eq. (2) (cf. Ref.^'®) 



(3) 



Depending on the ratio between the temperature of the 
system, T, the frequency under discussion, uj, and the 
applied voltage, one may distinguish among differ- 
ent regimes. In the limit {eV ^ maxjfcr, w}) the noise, 
eq.(3), is essentially non-equilibrium (shot noise), and is 
proportional to the absolute value of the d.c. electric 
current. For eV <^ maxj/cT, w} 52 is mostly thermal. 



or Nyquist noise. These observations are not system de- 
pendent, but rather follow from general considerations, 
namely the properties of various operators under time 
reversal transformation. 

Since the current operator changes sign under time re- 
versal transformation, any even-order correlation func- 
tion of the current fluctuations (e.g. ^2) taken at zero 
frequency is invariant under this operation. Assuming 
that current correlators are functions of the average cur- 
rent, /, it follows that even-order correlation functions 
depend only on the absolute value of the electric current 
(and are independent of the direction of the current). In 
the Ohmic regime this means that even-order current cor- 
relation functions (at zero frequency) are even functions 
of the applied voltage. Evidently, this general observa- 
tion agrees with the result eq. (3) for the second-order 
correlation function ^2 in the diffusive junction. 

By contrast, odd-order current correlation functions 
change their sign under time reversal transformation. In 
other words, such correlation functions depend on the 
direction of the current, and not only on its absolute 
value. Therefore in the Ohmic regime, odd-order corre- 
lation functions of current are odd-order functions of the 
applied voltage. This condition automatically guarantees 
that odd-order correlation functions vanish at thermal 
equilibrium. 

The present analysis focuses on the simplest, yet non- 
trivial, example of manifestly out-of-equilibrium correla- 
tor, namely the third order current correlation function 
S3. This is the lowest order correlation function which 
is dominated by non-equilibrium fluctuations even in the 
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regime where the applied bias is small compared with 
temperature. There are two situations where the study 
of shot noise under such conditions is called for. Firstly, 
when the lowering of the temperature of the electron is 
not facilitated. Secondly, and more interestingly, when 
one is interested in exploring the correlations in the sys- 
tem above a certain critical (or a characteristic) tempera- 
ture. On one hand this requires to keep the temperature 
relatively high. On the other hand, to study low energy 
features, one is restricted to low values of the voltage. 

In this work we argue and subsequently show that for 
the third order (and in fact all odd order) current cu- 
mulant, out-of-equilibrium fluctuations are not masked 
by large equilibrium noise: ^3 remains linear in volt- 
age (and temperature independent) at high temperatures. 
The first part of our analysis (Sectionll) focuses on the 
temperature dependence of . in the limit where the in- 
elastic length is much larger than the system's size. We 
find that 1S3 is proportional to the the d.c. current at 
all temperatures, and calculate how the proportionality 
coeflacient interpolates between the two asymptotic val- 
ues, 1/15 and 1/3, for T <^ eV and T > eV" respectively. 
In the second part of this paper (Section III) we study 
the effects of inelastic electron-electron scattering. We 
specifically consider the limit of hn/L <^ 1, where 1-^ is 
the inelastic electron-electron mean free path and L is 
the length of the conductor. We find that as compared 
with the elastic case 53 is weakly suppressed in the high 
temperature limit (by ~ 15%)but is highly enhanced in 
the low temperature limit (cf. Table I). Similarly to the 
elastic case ^3 is not masked by thermal fluctuations in 
the high temperature limit. 

Our analysis employs the recently developed non-linear 
(T- model -Kcldysh technique*~^°. This approach has been 
applied to study non-equilibrium noise in the presence of 
disorder and electron-electron interactions^*^. The main 



steps to be followed in the analysis below include finding 
the saddle point of the effective action (eqs. 18 and 43 
for the elastic and the inelastic cases respectively) and 
expanding in soft modes around them (eqs. 22 and 80). 
The saddle point (or the approximate one in the inelastic 
case) are connected with the single-particle distribution 
function, while the soft modes describe the dynamics of 
the density fluctuations. Both are significantly modified 
by strong inelastic scattering. 

There are of course other methods that can be em- 
ployed to study non-equilibrium fluctuations. One ap- 
pealing candidate would be the Shulman-Kogan^^ ver- 
sion of the Boltzmann-Langevin scheme. This method 
work successfully when ^2 is studied. In Section IV we 
present a short discussion which shows, arguably in an 
unexpected way, that the traditional kinetic equation ap- 
proach breaks down when higher order cumulants are ad- 
dressed. Section V includes a brief discussion. 



II. THE ELASTIC CASE 

As has been shown in Ref.^^, by modeling the measure- 
ment procedure as coupling of the system to a "quantum 
galvanometer" , the emerging mathematical object to be 
studied is 

Ss{x,t-,x\t'-,x'' ,t'') = {Tj2{x,t)l2{x' ,t')l2{x'' ,t'')) . (4) 

Here a; is a coordinate measured along a quasi one- 
dimensional wire (0 < a; < L) of cross-section A: Tc 
is the time ordering operator along the Keldysh contour 
and I2 is a quantum component of the current. Since the 
process we consider is stationary, the correlation function 
depends only on differences of time, — t'; t' — t"). In 
Fourier space it can be represented as 

(5) 



S3(u;i,a;2)= d{t - t')d{t' - t")e 



-ILJi {t — t')—iuj2{t'—t 



"^S3{t-t',t'-t") . 



(Here we have assumed that the frequencies cji , cl'2 arc 
small compared with the inverse diffusion time along 
the wire, such that the current fluctuations are inde- 
pendent of the spatial coordinate). We next evaluate 
the expression, eq. (4), for a particular system, in the 
hope that the qualitative properties we are after are not 
strongly system dependent. In the present section we 
consider non-interacting electrons in the presence of a 
short-range, delta correlated and weak disorder poten- 
tial (e/r ^ h, where r is the elastic mean free time 
and e/ is the Fermi energy). To calculate S3 we employ 
the (7-model formalism, recently put forward for deal- 
ing with non-equilibrium diffusive systems (for details 
see Ref.^°). As was stressed in our previous work this 
approach is a generalization of the Boltzmann-Langevin 
kinetic approach^ ^. It is comparable to the latter as long 
as the kinetic approach is applicable to diffusive systems. 



A direct employment of the kinetic approach, assuming 
the random force term to be short-ranged correlated in 
space, does not give correctly the the higher-than-two 
cumulants of the current. The reasons for that are dis- 
cussed below (cf. Section IV). 

The Hamiltonian of our system of non-interacting elec- 
trons in disordered systems is: 



Ho= [ 



dr^{r) 



Volume 



-|-(V-a)' + [/dis 
2m 



*(r). (6) 



Here ca/e is a vector potential. The disorder potential is 

(5-correlated: 



(«7dis(r) f/dis(r')) 



5{v-v'), 



(7) 



where v is the density of states at the Fermi energy. 
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Following the procedure outlined in Ref."'^'', employ- 
ing the diffusive approximation and focusing on the out- 
of-equilibrium system, one can write a generating func- 
tional, expressed as a path integral over a bosonic matrix 
field Q 

Z[a] = j PQexp(i5[Q,a]) . (8) 
Here the integration is performed over the manifold 

Q{x, t, ti)Q{x, ti , t')dti = d{t - t'), (9) 



the effective action is given by 
iS[Q, a] = -^Tr{D (VQ + i[a„7", Q]f 
and 
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, 72= (J 





4ieQ} , (10) 



(11) 



Tr represents summation over all spatio-temporal and 
Keldysh components. Here ai and a2 are the Keldysh 
rotated classical and quantum components of a. Here- 
after we focus our attention on ai, a2, the components 
in the direction along the wire. The third order current 
correlator may now be expressed as functional differen- 
tiation of the generating functional Z[a] with respect to 
a2 



83(^1—^2,^2 — ^3) = 



5^Z[a] 



8 &l2{xi,ti)Sa2{X2,t2)Sa2{X3,t3) 



■ (12) 



Performing this functional differentiation one obtains the 
following result 



83(^1- ^25*2- ^3) = 



^M{xi,ti)i^{X2,Xs,t2,t3) + 

(.Ti,ti ^a;3,i3) -I- {xi,ti ^X2,t2) + 



-M(a;i,ii)M(a;2,t2)M(a;3,i3) 



(13) 



P{x, x', t, t') = Tr''I^Q^,t,t'l2Qx',t',tl2 - St,tni |'^c.,x', (14) 
M{x,t) = Tr4 / dtJ[Q,,tM-M Qx,tut]l2]. (15) 



We employ the notations Q{x,t,t') = Q{x,t,t'); Tr^ is 
the trace taken with respect to the Keldysh indices; () 
denotes a quantum-mechanical expectation value. The 
matrix Q can be parameterized as 



Q = Aexp{W) , 

where 

AW + WA = , 
and A is the saddle point of the action (10) 



A(x,e) = 



_fl 2F{x,e) 







-1 



(16) 



(17) 



(18) 



The function F is related to the single particle distribu- 
tion function / through 

F{x,e) = l-2f{x,e) . (19) 

The matrix Wx,e.e', in turn, is parameterized as follows: 



(20) 



Here we have defined 

{w{x,ei,e2)w{x' ,e3,ei)) = 2(27r)^(5(ei - e4)5{t2 - €3)0(0;, a;', ei - €2) 



It is convenient to introduce the diffusion propagator 

{-iuj + £>V^)D(a;, x'uj) = -^6{x - x') . (21) 

nni' 

The absence of diffusive motion in clean metallic leads 
implies that the diffusion propagator must vanish at the 
end points of the constriction. In addition, there is no 
current flowing in the transversal direction (hard wall 
boundary conditions). It follows that the component of 
the gradient of the diffusion propagator in that direction 
(calculated at the hard wall edges) must vanish as well. 
The correlation functions of the fields w, w are then given 
by: 



{w{x,ei,€2)w{x',e3,e4)) = -g{2TTfS{ei - €4)6(62 - €3) j dxiB^^^^^^^^^^VF^^^^^VF^^^^^D^^^^^^^^^, 



{'w{x, ei, e2)w{x', £3, £4)) = 



(22) 



It is the the third cumulant (the reduced correlation 
function) of the current, 5*3, that is of interest to us 

53 = S3 - 3/S2 + 2/3 . (23) 



To evaluate S3 one follows steps similar to those that 
led to the derivation of S2, see Ref.^°. If all relevant en- 
ergy scales in the problem are smaller than the transver- 
sal Thouless energy {Eth = D/L\, where Lt is a width 
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of a wire), the wire is effectively quasi-one dimensional. 
In that case only the lowest transversal mode of the dif- 
fusive propagator can be taken into account, which yields 



Xi - X2\ - Xi - X2 + 



2X1X2 



• (24) 



Here g = hvD. The electron distribution function in this 
system is equal to 



and vanishes at equilibrium. The zero temperature re- 
sult (high voltage limit) has already been obtained by 
means of the scattering states approach for single-channel 
systems^^, and later generalized by means of Random 
Matrix Theory (RMT) to multi-channel systems (chaotic 
and diffusive) "'^^ . In our derivation we do not assume the 
applicability of RMT. Our result covers the whole tem- 
perature range. We see that at low temperatures the 
third order cumulant is linear in the voltage 



F(x,e) 



-F 



eV 
~2 



1 



eV 
~2 



(25) 



The quantities F and D determine the correlation func- 
tions, eq. (22). We can now begin to evaluate S3, (c.f. 
eq. (13)), performing a perturbative expansion in the 
fluctuations around the saddle point solution, eq. (18). 
After some algebra we find that in the zero frequency 
limit the third order cumulant is given by 



S3{U)1 = 0,W2 = 0) 



f 

J —I 



f 

Jo 



dx\dx2 



rfeF(e,Xi)D[0,Xi,X2]V(F'(e,X2)) . 



(26) 



Integrating over energies and coordinates we obtain 
53(a;i=0,W2 = 0) = e2/y(p) , 



y{p) 



6(-l -I- e^P) + 26e^P + e^P)p 
15p(-l -f e2p)2 



(27) 



where p = eV/2T. The function y is depicted in Fig.l, 
where it is plotted on a logarithmic scale. 
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(28) 



At high temperatures the electrons in the reservoirs are 
not anymore in the ground state, so the correlations are 
partially washed out by thermal fluctuations. One may 
then expand y{p), eq. (27), in a scries of eV/2T. The 
leading term in this high temperature expansion is linear 
in the voltage 



(29) 



Note that although thermal fluctuations enhance the 
noise (compared with the zero temperature limit), 
eqs.(28) and (29) differ by only a numerical coefficient. 
For the symmetry reasons alluded to in the introduction 
(the vanishing of iSa, iSs, etc. at equilibrium), it is tempt- 
ing to conjecture that all odd order current cumulants in- 
terpolate between two constant values as the ratio eV/ kT 
is varied. The experimental study of these cumulants 
provides one with a direct probe of non-equilibrium be- 
havior, not masked by equilibrium thermal fluctuations. 
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III. INELASTIC CASE 
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FIG. 1. The scaling function y plotted on a logarithmic 
scale, of. eq. (27) 

Let us now discuss the main features of the func- 
tion 5*3. In agreement with symmetry requirements (S3 
is an odd function of the voltage (the even correlator 
<S2 is proportional to the absolute value of voltage). 



In our analysis so far we have completely ignored in- 
elastic collisions among the electrons. This procedure 
is well justified provided that the inelastic length greatly 
exceeds the system's size. However, if this is not the case, 
different analysis is called for. To understand why inelas- 
tic collisions do matter for current fluctuations, we would 
like to recall the analysis of ^2 for a similar problem. The 
latter function is fully determined by the effective elec- 
tron temperature. Collisions among electrons, which are 
subject to an external bias, increase the temperature of 
those electrons. This, in turn, leads to the enhancement 
of <S2, cf. Refs.^^'^^. In the limit of short inelastic length 

kn^L, (30) 

the zero frequency and zero temperature noise is 

52(0) = y^el . (31) 

By comparing with the a;,T = limit of eq. (3) we see 
that inelastic collisions enhance the current fluctuations. 
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In the present section we consider the effect of inelastic 
electron collisions on 53. 

In the present analysis wc assume that the electron- 
phonon collision length is large, ^e-ph ^ L, hence 
electron-phonon scattering may be neglected. The 
Hamiltonian we are concerned with is 



H = Hq+ Hint- 



(32) 



The Coulomb interaction among the electrons is de- 
scribed by 

Hint = \j dYdY'^{v)^{v')Vo{v - r')*(r)*(r') , (33) 
where 



Vb(r-r') 



(34) 



Wc need to deal with the effect of electron-electron inter- 
actions in the presence of disorder and away from equi- 
librium. Following Ref.^ one may introduce an auxiliary 
bosonic field 



$ = 



(35) 



which decouples the interaction in the particle-hole chan- 
nel. Now the partition function (eq. (8)) is a functional 
integral over both the bosonic fields Q and (j) , 



The action is 



Q2 = l 



'DQ'D(j)ex.p{iStotai) ■ 



iStotai=iSm+iS[^,Q], 
iS[^] = iTr{$'^yo"S^*} , 

iS[^,Q] = -^Tr{Q2}+Trln 

4r 



iQ 
2t 



i7 



(36) 



(37a) 
(37b) 

(37c) 



It is convenient to perform a "gauge transformation" ® to 
a new field Q 

Qt,t'{x) = exp {ika{x,t)Y') Qt,t'{x) cxp {-ika{x,t')-f") . 

(38) 

Introducing the long derivative 

d^Q = VQ + i[Vkaj",Q], (39) 

one may write the gradient expansion of eq. (37c) as 

iS[Q, = iuTi{{^ - iujKfj2{^ + iojK)} - (40) 
DTr{d^Qy + 4iTr{(e + {4,^ + iLok„)-f'')Q} . 



TTV 

T 



At this point the vector = (fcijfe) that determines 
the transformation (38) is arbitrary. The saddle point 



equation for Q of the action (40) is given by the follow- 
ing equation 

Dd^iQd^Q) + i[ie + (</.„ + iLu%,h''),Q] = . (41) 
Let us now choose the parameterization 

Q = Aexp(H^), (42) 
where W represents fluctuation around the saddle-point 
'1 2F[<j)]{x,ey 



A{x,e) 







(43) 



Eq. (43) implies that the solution of the saddle point, 
equation (41), determines F as a functional of </>. We 
do not know, though, how to solve it. Instead we aver- 
age over cf) on the level of eq.(41). The solution of this 
averaged equation, denoted by F, is determined by: 



DV'^F{e) = r^lF} 



(44) 



where the r.h.s. is given by 

7^«{^} = D J ^[{Vk\uj)Wk\-uj)){F{e)-F{e-uj)) + 

{{Vk^ {uj)Vk''{-uj)) - (^2(^)vfci (-^)))(F{e)F{e- oj) - 1)] . (45) 

We stress that our averaging procedure is not related to 
the genuine interacting saddle point. We will be expand- 
ing about another point (A) at which the fields W and 
(f), (Vfc) are nearly decoupled; W is parameterized as in 
eq.(20). Taking variation of the action with respect to w, 
w, we obtain the following gauge, determining k[(j)]: 



DV^k2 -4)2- iojk2 = 

DV'^ki + (j>i+ iLoki = 2B[u!, x] V^A;2 



(46) 



where 



B[a;, x] 



F{e,x)F{e-uj,x)]. (47) 



Though we have failed to find the true saddle point the 

linear part of the action expanded around (A) is zero. 
It is remarkable to notice that under conditions (46) 
eq.(44) becomes a quantum kinetic equation^^ with the 
collision integral being The above procedure 

consists of an expansion of the action, eq.(40), aroimd 
its saddle point. The expansion presumes slow fluctu- 
ations in space and time (though the results obtained, 
expressing ^3 in terms of the distribution function and 
the effective propagator, eq.(79), may be valid imder less 
restrictive conditions). In Fourier space (with respect 
to the spatial coordinate and the sum-rather than the 
difference- of the time coordinates) this implies ql <^ 1 
and WT <C 1. Hence this procedure implies that the 
collision integral, 7^^, appearing in the quantum kinetic 
equation, eq.(44), accounts only for electron-electron col- 
lisions with small momentum transfer. We therefore need 
to ascertain that these are indeed the small momentum 
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transfer processes which make the dominant contribu- 
tion to the electron-electron collision rate. To estimate 
the relative importance of the small and large momen- 
tum processes we compare their respective contributions 
to the energy relaxation time, Tee- Electron-electron scat- 
tering with small momentum transfer leads to result^^'^^: 



e < 



(50) 



1/Tee = 



,d/2 



(48) 



while the contribution coming from large momentum 
transfer is 



V-^ee = ^ » 1 . 



(49) 



the contribution to the inelastic mean free time coming 
from scattering with small momentum transfer is domi- 
nant; ignoring it would result in seriously overestimating 
T^^ time (and correspondingly the mean free path Zee)- 

Coming back to our calculations we note that the cor- 
relation function of current fluctuations is a gauged in- 
variant quantity (does not depend on the position of the 
Fermi level). This means that momenta q < ^/ujJD do 
not contribute to such a quant ity -'^ . In this case the 
Coulomb propagator is universal, i.e. does not depend 
on the electron charge. The fact that we address gauge 
invariant quantities allows us to represent the generating 
functional Z in terms of the fields Q and V/c (rather than 
Q and <p), as in Ref^. 



As is now clearly seen, for low enough energies 

{Z) = jvVK exp (-i!/£)Tr{Vii''^p-^V/r}) j VQ exp f V/C] 



\i=0 



Here we define 



-2iLuSx,x'B^{x) 



where the expansion 5 = 5° -|- 5^ -|- 5^, is in powers of \/K; the I — th power (Z = 0, 1, 2) is given by 



iS°[Q] = - ^ [DTr{VQy + Ai Tr{eQ} 
iS^lQ, VK] = -inu \DTv{Vkaj"QVQ} + Tr{((?i„ + iLokc,)^"^} 



Tr{VA:„7"gV/c;37^<9} - Tr{VA;«7"AVA;^7^A} 



Prom eq.(51) we obtain the gauge field correlation function 



{\/ka{x,U!)\/k0{x',-Ul)) = —Ya^/s{u!,X,x'), 



where 



Y{clJ, X, x') 



—2i'KVUJ J dx\T)[—Lo,x,xi\Q[LjO,x-i\'D[oj,x\,x'] D[— a;,.T,.r'] 
T)[ui,x,x'\ 



Using eqs. (54,55) we rewrite eq.(45) for the quasi-one-dimensional wire as: 
£)V^F(e) =/""{e,a;},where 



P^{e,x) = "-^ j dw[-2i(jjTTuD[x,xi,-uj]B[uj,xi]T)[xi,x,(jo]{F{e) - F{e- uj)) + 
{T>[x, x, w] - D[x, x, -L0]){1 - F{e)F{e + lo))]. 



(51) 
(52) 

(53a) 
(53b) 
(53c) 

(54) 

(55) 

(56) 

(57) 



The total number of particles and the total energy of 

the electrons are both preserved during electron-electron 
and elastic electron-impurity scattering. The collision in- 
tegral, eq.(57), satisfies then 



r"(e,a;)de = 0. 



/OO 
er^(e,x)de = 0. 
-OO 



(58) 
(59) 
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Wc now consider the limit ^cc ^ L. The sohition of 
eq.(44) assumes then the form of a quasi-equihbrium 
single-particle distribution function 



F{e, x) = tanh 



e — e(f){x) + /X 



(60) 



Here e is the total energy of the electron, and e(j> is the 
local electro-chemical potential, i.e. the sum of the lo- 
cal chemical potential /z and the electrostatic potential 
u {e(j){x) — iJ.{x) + e?i(.x)); p, is the value of the electro- 
chemical potential before the bias has been applied and 
T{x) is the effective local temperature of the electron gas. 

Since the size of a constriction is mu ch longer than 
Debye screening length A^^ = V Ane'^i' we may use a 
quasi-neutrality approximation, in this case to assume 
that a value of chemical potential /U is constant in the 
constriction. It still remains to find the effective elec- 
tronic temperature, T{x), and the electrostatic potential 
u{x). 

In order to find the electrostatic potential that enters 
eq.(60) wc employ eq. (58). To facilitate our calculations 
we further assume that conductance band is symmetric 
about the Fermi energy and that the spectral density of 
single-electron energy levels is constant. Integration over 
the energy, eq.(58) yields 



dl / F{e)de = 



= 0. 



(61) 



(62) 



Solving eq.(62) with the requirement that the electro- 
chemical potential at the edges of the sample differs by 
V, we find that the electro-chemical potential along the 
constriction (0 < x < L) is given by 



e(l>{x) = eV[^-^] +Jx. 



(63) 



We use the energy conservation property of the colli- 
sion integral, eq. (59). Multiplying eq. (56) by energy 
and integrating over it 



: j eF{e)de = 







(64) 



we obtain an equation for the electron temperature: 



dl(^{kT{x)r + \{euix)r]=0. 



(65) 



The boundary condition of eq.(65) is determined by the 
temperature of the electrons in the reservoirs. Combin- 
ing eqs.(64) and (63) we find the electron temperature in 
two opposite limits: 



T{x) 



^^x{L-x) eV-»T 
T eV 



(66) 



Eqs. ((60) , (63) and (66) determine the fimction F 
uniquely. We now replace the right-corner element of the 
matrix A (i.e. cf. eq.(43)) by its average value F. 

To calculate iS,-? under conditions of strong electron- 
electron scattering (eq. (30)) one needs to replace the 

operators I and M in eq. (13) by their gauged values 

^3(^1 — t2, t2 — t3) = 

-M{xi,ti)l (x2,X3,t2,h) + 
{Xi,ti ^X3,t3) + (Xi,ti ^X2,t2) + 
^^^^A{xi,ti)A{x2,t2)Mx3,t3)\ 



(67) 



where the averaging is taken over the entire action S and 
the Gaussian weight function for \7 K, as in eq. (36). Here 
we define (cf. eqs.(14), (15) with eqs. (68), (69)) 

i {x,x',t,t') =Tr\Q^^t,t'l2Qx',t',t'y2 - St,tm\Sx,x' , (68) 



M{x,t) = Tri^J dh(^[Q^,t,t,;d^] QxM.t^l2j , (69) 

where the "long derivative", dx, is presented in eq.(39). 
In order to actually perform the functional integration 
over the matrix field Q we use the parameterization of 

eq. (42). The operators I ,M may be expanded over 
w, iu and VA; : 



M = Mq + M? + . . . + Mo + m| . . . 

= io + ii + . . . 



(70) 
(71) 



Here the upper index refers to the power of the Vfc 
fields: the lower refers to the power of u\ ill fields in the 
expansion. We need to find the Gaussian fluctuations 
around the saddle point of the action (53). Though we 
did not find the exact saddle point, the expansion of Q 
around A works satisfactorily. The coupling between the 
fields V/c and W which appears already in the Gaussian 
(quadratic) part is small, since it is proportional to the 
gradient of the distribution function: 



(72) 



Considered as a small perturbation, iSj does not affect 

the results. 

The more dramatic effect on the correlation function 
arises from the non-Gaussian part of the action, eqs. 
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(53b, 53c) (by this we mean non-Gaussian terms in ei- There is no hncar term in Vfc and the Gaussian part has 
ther or WK). After integrating over the interaction been separated out the following identity hold: 



an additional contribution to the Gaussian part (propor- 
tional to ww) of the action arises. To find the effective 
action iS^^ [W] we average over the interaction along the 
following lines: 



exp <^ iS^ + iS"' 



(74) 



Vk 



exp I + iS'^) + +iS'^-{ iS^ + iS^ 
exp \(iS^ + iS^ 



(where, again, Sq refers to the component of the action, 
eq.(53), that has zero power of the w,iju fields and one 
power of the Vfc field). In addition, due to the choice 
of the gauge, eq.(46). and the condition (Zee ^ L), the 
averaging over Vfc does not generate terms linear in w, w 
in the effective action: 



VA; 



{iSf} 



Vk 



/d€. 



(75) 



Vk 



+ iS'^-(iS^+iS'^j +-\iS^+iS^-(iS^+iS^ 



^ Combining eqs. (73,74 and 75) we find that the effective 
action acquires an additional contribution: 



:^expl^(iS^+iS^^^(j^iS^ + iS^^ "^-^(iS^+iS^^y (73) (^exp {iS^ + iS^)^ cjip (^{iS^) + ^{iSfiS^)^ (76) 

The general form of the effective action is rather complicated, however for the low frequency noise only diagonal 
part of the action matters: 



iSf[w,xS]= — TV<jiD,,, 



DV^+X^ 



(77) 



Here the operator 



(78) 



is a linearized collision integral, i.e. a variation of the collision integral (57) with respect to the distribution function. 
Inspecting eq.(67) we find (details are outlined in Appendix A) 



S-i{uj\ = 0,a;2 = 0) = 



Se^Ang 
nL3 



2 fL 



I 

Jo 



dx\dx2 



/oc 
(iei(ie2F(ei,a;i)I?[xi,ei;x2,e2]V (-^^(£2, 3^2)) • 
-00 



The diffusion propagator in the presence of strong inelastic electron-electron scattering is given by^'^ 



2?[a;i,ei;x2,e2] = ( ^/o 



D[a;i,a;2] ( -1 



3 ei - n{xi) 02 - fi jxi) 
7r2 T{xi) T{xi) 



(79) 



(80) 



where 



k{x) = 



1 



1 + exp(a;) 

Evaluating .Ss explicitly we find that the third order current cumulant is 



S3{lOi =0,U)2 = 0) 



36e^Ag'^eV 



dxidx2D[xi, X2] 



T{x 



T{X2) 



-I- (a:i - X2)-^ 



d 



T{X2) dx: 



(81) 



At high temperatures (cf. eq.66) one obtains 
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53(a;i =0,a;2 = 0) = ^e^/, 



(82) 



while at low temperatures 

53 (wi = 0, W2 = 0) = (^-^ - e^I. (83) 

Our calculation was performed for a simple rectangular constriction. However, our results hold for any shape of 
the constriction, provided it is quasi-one dimensional (we have considered a single transversal mode only). 

Finally we discuss the role of an applied magnetic field. Consider a two-terminal geometry and a non-interacting 
electron gas. For such a system any current correlation function can be expressed through the channel transmission 
probabilities {T^}. Onsagcr-typc relations dictate that the transmission probabilities arc invariant imdcr the reversal 
of the sign of the magnetic flux ("phase locking"), i.e., Ti(<^) = Ti{—(j)). This implies iS3(0) = <S3(— </>). Let us take 
the zero bias voltage limit. Because S3 is an even function of the magnetic flux it cannot change its sign under 
time-reversal operation. On the other hand, by its very definition, it must reverse its sign under the time reversal 
operation. It therefore must vanish. We conclude that the magnetic field alone, while breaking the time reversal 
symmetry, cannot lead to finite-value odd current correlators. 



IV. WHY THE STANDARD KINETIC EQUATION APPROACH FAILS 

The phase coherence is not significant for the current fluctuations in the weakly disordered metals (the statement 
remain true for any system for wliic;li electron dynamics is classical) . This suggest that the kinetic theory of fluctuations 
may be suitable description of the problem. Here we state, that even though the dephasing length may be short 
{Irmph ^ L) the is no purely classical route to flnd high order current correlation functions. 

One possible way to represent the kinetic theory of fluctuations is through Boltzmann-Langevin equation^ ^. The 
deviation 6f of the exact distribution function from its coarse-grained) value (/) satisfles a stochastic equation with 
the random additive noise {SJ): 

+ V • v) <5/(p, r, t) = Col{/} + 6r-\p, r, t). (84) 

The l.h.s. of eq. (84) is the conventional drift term of the kinetic equation; Col is a collision integral. The random 
flux term SJ which is added to the standard Boltzmann equation, induces fluctuations of the distribution function /. 
The correlation of the random fluxes can be related in the universal manner to the coarse-grained particle flux in the 
phase space (which is specific for each given problem). However the very applicability of the kinetic equation as well 
as the applicability of kinetic theory of fluctuation arc controlled by the same parameter (in our case epT). However, 
the normalized ratio between high and low order current cumulants is small by exactly the same parameter. 

To see, where the kinetic (Kogan-Shulman) procedure fails we try to go along the standard way (for the simplicity 
for non-interacting electrons). In the diffusive approximation one may keep only the lowest harmonics: 

Sf{p, r, t) = Sfo ip, v,t) + n- 5^ {p, r, t) , (85) 



W«^*(p, r, t) = SJoip, r,t) + n-SJi {p, r, t), (86) 

where p ~ n\pp\ (here n is unit vector in the direction of the momentum). Here fi and Ji are vectors in the direction 
of the average current (the x-axis). From eqs. (85,86) one flnds 

6h = -vftVSMp, r) - tSMp, r, t). (87) 

In the low frequency limit particle conservation guarantees that current fluctuations through any given cross section 
are the same. This allows us to compute the current fluctuations integrated over the entire sample length, 

Sm = JdrJ de63^,{e,r,t). (88) 

In order to calculated the third order current cumulant 53(0,0) it is crucial to know the zero frequency limit of the 
third order correlation function of the random flux. 
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7^3(l"l,r2,r3, 61,62,63)= J d{ti - t2)d{ti - t^) 

((5Jia;(ei, ri, ti)SJi^{e2, r2, i2)'5Jix(e3, r3, t3))- (89) 

If one would go along the standard assumption (of Poissonian, and ^-correlated in space random flux), one will get 
a wrong answer (smaller as (l/L)^ than the correct one). Instead one can show, that by choosing the correlation 
function to be long correlated: 

27 TT 

^3(ri,r2,r3,ei,e2,e3) = — - — r^(5(ei -e2)(5(ei -e3)5(ri -r2) 

S{yr - y3)6{zi - zs){\/F\euXi))D{xuX3)F{ei, xs) (90) 

eq. (26) (and similarly eq. 79)) are reproduced. Here the function F{F) is related to the singe-particle distribution 
function (cf. eq. (19)). Eq.(90) shows that the correct high-order correlation functions of the random fluxes are 
non-local in space. This is clearly in contradiction with the assumption made within the Kogan-Shulman formalism. 

In other words, one usually employs the Boltzmann-Langevin-Kogan-Shulman approach approximating the actual 
probability distribution function of the random current flux by a Gaussian distribution. Seemingly, this is justified 
because of the large dimensionless conductance <? » 1. In the scattering states picture, this condition is equivalent to 
a large number of transmission channels. We know^^ that each channel provides its own independent contribution to 
the current noise. According to the central limit theorem, a process that consists of a large number of independent 
random contributions is characterized by a (nearly) Gaussian distribution function, independently of the distribution 
function of the individual contributions involved. However, because the value of the dimensionless conductance is 
indeed large, yet finite, the distribution function of the random current fiux deviates from Gaussian. Such deviations 
are beyond the validity of the Boltzmann-Langevin approach. S3 is the lowest order current cumulant that probes 
these deviations. 

This is not surprising. Indeed, since the distribution function / is a macroscopic quantity, (at equilibrium) it must 
satisfy the Onsager regression hypothesis^^ . The kinetic theory of fluctuations can be viewed as an extension of 

the this hypothesis for a non-equilibrium situation^^. Since the Onsager hypothesis is restricted to pair correlation 
functions only, it is natural to expect that the kinetic theory has the same limitations. 



V. DISCUSSION 



We have noted that due to general symmetry reasons odd order cumulants of the current must be odd functions of 
the applied voltage and must therefore vanish at equilibrium. Provided these cumulants arc analytic functions of V, 
they must increase linearly with voltage bias when the latter is small, they increase linearly with the latter. Thus, 
contrary to even order cumulants, they are not masked by a background thermal noise. Even more remarkably, at 
high temperature, the third order current cumulant approaches a constant value. This is in contradistinction with 
the second order current correlator which diverges with the temperature. This is the reason why quantities such as 
Ss are suitable for the study of non-equilibrium current noise even at relatively high temperatures. 

Our results for .Ss, pertaining to both the elastic and the inelastic cases (the high and low temperature limits), are 
summarized in Table I: 





T< eV 


T^eV 


kn > L 


1/15 e^7 


1/3 e^I 


hn < L 


0.248 e^I 


0.304 e^I 



TABLE I. Limiting values of S3. 



In the elastic limit one can regard the electrons as non-interacting; the problem then is a natural extension of 

the ballistic multichannel setup. Indeed, at low frequencies one may ignore the dynamics of the electrons within the 
constriction, which allows us to employ the Landauer scattering states approach^*. The latter consists of describing 
the constriction by a quantum transmission matrix, expressing the currents and their correlation functions through 
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the transmission coefficients. (It has been noted, though, that coherency is not essential for deriving the leading order 
contributions to shot- noise in disordered systems, cf. Refs.^'^'^'*). 

Since odd order current correlation functions vanish at equilibrium, they must be proportional to the difference 
between the distribution functions on left and right respectively. But such contributions do not diverge with increasing 
temperature, which is a formal reason for the saturation of the odd order correlation functions at high temperatures. 

For disordered systems the transmission coefficients arc random variables and their statistical properties are known 
from Random Matrix Theory (RMT)^^'^^. Employing RMT one may average the scattering states result over the 
disorder. This reproduces the results we have obtained in the limit hn L, and in particular eq.(27). 

Finally, we briefly comment on the observability of high order current correlation functions. Consider an ametcr 
that detects the net charge transmitted through a given cross-section over a time interval t. Under the condition 
that r is sufficiently long (n = (/)r/e » 1, where n is the average number of electrons passing through a given 
cross-section within the time interval r) , one may extract the low frequency current fluctuations from the statistics of 
the transmitted charge (known as a counting statistics). Indeed, to describe low frequency current fluctuations one 
has to require h/r to be smaller than the relevant energy scale. For temperature T ~ lOOK we require t > 10~^p sec, 
the latter condition satisfied for any practical measurements. The more realistic restriction comes from the available 
electronics, and we estimate r ~ Insec. For (/) ^ 100 nA we find n ~ 100. 



(((n - n)2))3/2 ni/2 lo' 

As we see, in principle a measurement of S3 is possible. One should note that in this consideration the nongeneric 
parasitic noise (such as an amplifiers 1// noises and etc.) had been ignored. 

Zero temperature current noise in quantum coherent junctions (i.e. the elastic case)"'^'* is a manifestation of a 
stochastic process that follows binomial distribution. The full counting statistics in the elastic limit has recently 
been found"'^^'^'''^°. It remains an open problem to find the counting statistics in the strongly inelastic limit. Our 
analysis reveals that the asymmetry of the probability distribution function (a measure of which is S3) is modified 
in a non- trivial way (i.e., it is either enhanced or suppressed by the electron-electron interaction, depending on the 
dimensionless parameter eV/T). 

Upon completion of this paper, we have learned of a related work by Levitov and Reznikov (cond-mat/0111057) 
addressing similar questions in quantum coherent conductors. We acknowledge discussions with M. Reznikov, A. 
Kamcnev, A. Mirlin and Y. Levinson. This research was supported by the U.S. -Israel Binational Science Foundation 
(BSF), by the Minerva Foundation, by the Israel Science Foundation of the Israel Academy of Arts and Sciences and 
by the German-Israel Foundation (GIF). 



APPENDIX A 



In this appendix we present the technical steps involved in evaluation of the 53 in the <JC L limit. Starting with 
eq.(67) one notes that the part which contribute to the correlation function at question is given by: 

Ssitx - t2M - h) = / Ml (a;i , 1 1 )Ii (a;2 , X3 , i2 , is ) + ^-Mi (a;i , i 1 )M2 (x2 , ^2 )Mi (xa , ) -h 

Mo(xi,ii)if (x2,X3,t2,i3) +M°(xi,ii)Mj(x2,i2)Mo(xa,i3 + {xi,ti ^xa.ia) + {xi,ti ^X2,h)^ . (Al) 

The first two terms entering eq.(Al) are similar to those we dealt with in the non-interacting case (cf. (13)). The 
third term in eq.(Al) should be kept because out of equilibrium the fields Vfc and w are no longer decoupled (eq. 72). 
However, since the coupling between those two fields is proportional to VF we may expand the action in the latter: 

Ssih - t2, t2 - ia) = (M^{Xi,ti)li {X2,X3, ^2, +^^M^{xi,h)M2{X2,t2)M^{x3, t3) + 

Ml{xi,ti)lf {X2,X3,t2,t3)iSl +Ml{xi,ti)l^^^^^ + {Xi,ti ^X3,t3) + ixi,ti ^X2,t2)) ■ (A2) 

/ 

Here ()o implies that the expectation value is calculated employing the correlators eqs.(54) and (22). The values of 

-D 

I and M are then given by; 
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i°{x,t2,t3)=-2 j d[e]i 



,i(t2(e3-e6)+t3(e5-«4) 



e5,€6,x"e3,e4j ) 



Mi(x,t) = / d[e]e'('^-'^)'[Vfcu3_,,,^u;^,,,,e3-V/ci,,_,3,^«;^,,3,,J. 



(A3) 
(A4) 
(A5) 

(A6) 



To evaluate >S3(ti;i = 0,a;2 = 0) (eq (Al)) we will use 
the fact that the current fluctuations are independent of 
the choice of the cross-section; wc therefore may inte- 
grate eq.(Al) over the entire volume of the sample. We 
find that in the limit Zin 4C i (Ss has the structure (cf. 
eq. (26)) 

Ze^Ang"^ f°° 
SsiiVi = 0,u}2 = 0) = — — / dxidx2 / rfeirfe2 

F{e,x^){w^,,,^,,w^,,,^,,)V{F\e,X2)) . (A7) 

Here the expectation value is calculated accordingly to 
the effective action (77). Performing the Gaussian inte- 
gration we find the correlation function: 

(tt^x,ei,eiW'x,e2,e2) = 2D(a;i, Ci; X2, £2) • (A8) 

Here we define the zero frequency a propagator: 



-DV^ + I-'\D{xuei;x2,e2) 



^(£1,62) 



5{xx - X2). (A9) 



The propagator!) is nothing else, but the kernel (taken at 
zero frequency) of the kinetic equation (in the diffusion 
approximation). In the limit when electron scattering 
rate is small, the inelastic diffusion propagator becomes 
an ordinary diffusion. In the opposite limit, one can in- 
voke the perturbation theory in L/P" to find its explicit 
form. Together with eq.(A7) its leads us to eq.(79). 
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